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ABSTRACT 



Context. A number of accreting X-ray pulsars experience spectral changes, both on the long time scales and on the time scales of 
the neutron star spin period. The sources seem to form two distinct groups which dilfer by the type of the spectral variations with 
flux. Such a bimodality probably reflects two dilferent regimes of accretion that may result in a particular pulsar depending on its 
luminosity - so-called sub- and super-critical regimes. 

Aims. We investigated the behavior of the spectral parameters of the Be/X-ray binary system A 0535 +26, as a function of flux and 
pulse phase. 

Methods. We used the data collected with INTEGRAL and RXTE during the April 2010 outburst of the source. We analyzed the 
phase-averaged and phase-resolved spectra and performed pulse-to-pulse spectral analysis of the pulsar 

Results. Our analysis reveals variability in the continuum parameters of the source's pulse-averaged spectrum with flux. The pulse- 
averaged cyclotron line energy does not change with the source luminosity during the outburst, which is consistent with previous 
studies. Our pulse-phase resolved and pulse-to-pulse analyses reveal, however, indications for a positive correlation of the cyclotron 
line energy with flux, as well as a flux-dependence of the continuum parameters. Based on the observed behavior, we argue that 
A 0535+26 operates at the sub-critical accretion regime. 

Key words. X-rays: binaries - Stars: neutron - Accretion, accretion disks - pulsars: individual A 0535+26 
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1. Introduction 

1 The X-ray transient A 0535 +260 was discover ed by the satellite 
• Ariel V during a giant outburst in April 1975 (iRosenberp et al.l 
' Il975h . The source belongs to the class of Be/X-ray binaries. It 
consists of an accreting neutron star orbiting the 09.7IIIe op- 
tical star HDE 245770 (Li et al. 1979; GiangrandLe et al. 198§. 
The system regular ly undergoes X-ray outbursts separated by 
quiescence phases jMotch et al.l 1199 lb . In Be/X-ray binaries, 
outbursts are often divided into normal and giant ones. Normal 
outbursts repeat periodically at times near the periastron passage 
and are thought to originate from an increase of the accretion 
\ rate (Lmax/^'min ~ 100) due to a closer positio n of the neutron 
i star to its Be star companion (IStella et al.lll986l) . Giant outbursts 

■ show a larger increase in lumin osity (Lmax/^m in ~ 100-1000) 
and last several tens of days (St ella et al.l ll986'). Often there is, 

■ however, no strict distinction between the two types of outbursts. 
Thus, some outbursts in A 0535+26 can not be clearly assigned 
to one of those two groups. Besides the outburst in 1975, 
A 0535+26 w ent through gi ant outbursts in 1980 (Naaase et 
[19821) . 1983 (Sembay etal.'ll990|), 1989 (Makino et al. 098^ 
1994 Jpingeretal.. .19961 2005 dTueller et aL 2005,1 2009 
(iKrimm et alJ l2009t ICaballero et all l2009l l2011ch and 



201 1 ("Camero-Arran z et al. I 120111: lTcheminetal.ri2011t 
ICabal lero et al. 2011a). 

Th e distance of the system has been measured to be around 
2kpc dSteele et al.lll998l) . Recent measurements of the orbital 



' A 0535+26 is referred to as 1 A 0535+262 in SIMBAD 
Astronomical Database 



period give ~ 1 1 1 . 1 days ("Finger et al." 2006h . The pulse period is 
around 103.39 s (Caballero et al..2007D . 

The phase-averaged X-ray spectrum of A 0535+26 can be 
described by a hard power law (r = 0.8- 1.1 for 2-18keV 
Ariel-V/SSI observations. Rick etts et al.lll975h modified above 
~ 15-20 ke V by an exponential ro ll-off (TTM/HEXE spectrum in 
2-156 keV, iKendziorra et al.llT994 ). Alternatively. lAcciari et al.l 
(1201 ih used the thermal Comptonization model to fit the low- 
energetic part of the X-ray continuum. Additionally, the spec- 
trum shows absorption lines, interpreted as cyclotron resonant 
scattering features (CRSFs). Such features are caused by res- 
onant scat tering of photons off the electrons in Land au levels 
(see, e.g.. iTruemper et anil978l; Usenberg et al.lll998b . The en- 
ergy of the fundamental line and the spacing between harmon- 
ics are directly proportional to the magnetic field strength. Two 
cyclotron absorption lines at about 50keV and lOOkeV were re- 
ported and interpreted as the fundamental and the first harmonic 
and implying a magne t ic field strength of a bout 4.3 • lO'^G 
( IKendziorra etalj|1994l: ICaballero et"al]|2007l) . While the pres- 
ence o f the fundamental line at about 50 keV was initially ques- 
tioned (iGrove et al .ll995HArava & Hardingll996l) . newer analy- 
ses clearly confirm the presence of both lines and their interpre- 
tation as a fundamental cyclotron feature and its first harmonic 
(ICaballero et al.ll2007l) . 

An evolution of the cyclotron line energy with the source lu- 
minosity has been reported for several accreting X-ray pulsars. 
A negative correlation of the line energy with luminos ity has 
been observ ed in V 0332+53 and 4U 01 1 5+63 (Tsygank ov et al.l 
2006, 200l l20Tol iMowlavi et al.ll2006l) . A positive correlation 
was detected for HerX-1 and GX 304-1 (IStaubert etal.ll2007l 
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iKlochkov etalll2012h . A similar behavior is seen in the pulse- 
to-pulse variations. For V 0332-1-53 and 4U0115H-63, a neg- 
ative correlation of cyclotron line energy with the amplitude 
of individual pulses and a softening of the spectral continuum 
with pulse amplitude was reported, whereas for HerX-1 and 
A0535H-26 a positive coiTelation of the line energy with pulse 
amplitude and a hardening of the spectrum with pulse amplitude 
was observed (Klochkov et al. 2011). Since A 0535-1-26 showed 
the spectral-flux dependence on the pulse-to-pulse timescales 
only, it is important to verify and study this behavior, using the 
data from a new outburst of the source. 

Pulse-phase resolved analysis shows that the spectrum of the 
main peak of the pulse profil e seems to be harder than the spec- 
trum of the secondary peak (iKendziorra et al.|[l994t ICaballerol 
l2009l) . ICabaiierol(l20OTh also found a noticeably lower cyclotron 
line energy during the main peak compared to the secondary one. 

In this work, we analyzed the April 2010 outburst of 
A 0535+26. On April 3 2010, a brightening of A 0535 +26 in 
the Swift/BAT light curve was reported (Caball ero et al.l2010l) . In 
the 20-1 00 ke V energy band of INTEGRAIJISGRI the flux of the 
source reached a value of about 0.92 Crab and the fundamental 
cycl otron line was detected at about 47.6 keV (ICaballero et alJ 
12010) . The outburst analyzed in this work took place after 
the giant outburst of December 2009, which reached a flux 
of about 5. 1 4 Crab in the 15-50keV Swift/BAT fight curve 
( ICaballero e ~al. 2010). 

Our analysis is focused on tfie spectral evolution of 
A 0535+26 as a function of luminosity. We analyzed pfiase- 
averaged, pulse-pfiase resolved, and pulse-amplitude resolved 
(pulse-to-pulse) spectra of tfie 2010 outburst, to explore spectral 
variations on sfiort (down to tfie pulsar's spin period) and long 
(several days) time scales. A description of tfie observations and 
data processing is given in Section |2] In Section [3] we describe 
tfie analysis of tfie pfiase-averaged, pfiase-resolved and pulse-to- 
pulse spectra. In Sections |4] and |5] we discuss and summarize 
our results. 



2. Observations and data processing 

For our analysis, we used data from tfie JEM-X and IBIS instru- 
ments of tfie INTEGRAL satellite (INTEmational Gamma-Ray 
Astropfiysics Laboratory, WjnkleoJLal. 2003) and from tfie PCA 
and HEXTE i nstruments onboard Rossi X-ray Timing Explorer 
(RXTE. Bradt et al.lll993l) . 

Tfie JEM-X monitor provides spectra and imag es witfi arc- 
minute resolution in tfie 3-35 keV energy band (iLund et al.l 
[2003) . Tfie fiigfi angular resolution gamma-ray im ager IBIS con- 
sists of two detector layers, ISGRI and PICsIT dUbertini et al.l 
l2003h . Tfie ISGRI layer, wfiicfi we used for our analysis, is sen- 
sitive at energies from 15keV to 1 MeV (Lebrunetal. 2003). 
Tfie PCA instrument onboard RXTE consists of a set of five 
proportional count ers and is sensitive to energies 2-60 keV 
dJafioda et al.ll2006l) . Tfie fiigfi-e nergy detector HE XTE operates 
in tfie energy range 15-250keV (iRotfiscfiild et al.l[T998ft . 

We used data covering tfie time from tfie maximum to tfie 
decay of tfie April 2010 outburst (2010-04-03 to 2010-04-17, 
Table [B- A Swift/BAT ligfit curve of A 0535+26 witfi tfie indi- 
cated observations is sfiown in Fig.[T] 

Tfie INTEGRAL data were processed using OSA 9.0. For tfie 
ISGRI data we performed an additional gain correction to ac- 
count for known problems in tfie 0SA9.Q energy calibration. 
We used tfie JEM-X data in tfie 5-30 keV range and tfie ISGRI 
data in tfie 20-145 keV range. We added systematic errors of 
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Fig. 1: Tfie 5w//f/BAr daily ligfit curve of tfie April 2010outburst 
of A 0535+26. Tfie INTEGRAL and RXTE observations wfiicfi 
we used for our analysis are indicated witfi dasfied (INTEGRAL) 
and solid lines (RXTE). 

Table 1 : Selected observations. 



Observation Exposure [ks] Flux in Crab [10- 100 keV] 



INTEGRAL0912 65.1 0.978 

1NTEGRAL0914 23.8 0.675 

INTEGRAL0915 28.6 0.368 

INTEGRAL0916 29.4 0.191 

INTEGRAL 0917 3L7 0.096 

RXTE 95347-02-01-00 3.3 0.908 

iUCre 95347-02-01-01 3.3 0.986 

RXTE 95347-02-01-02 3.0 0.715 



3% for JEM-X and 2% for ISGRI, wfiicfi we estimated by fitting 
a Crab spectrum accumulated witfi tfie respective instruments 
witfi a canonical power-law witfi a pfioton index of ~2.1. Tfie 
Crab spectra are taken during tfie same observations since botfi 
sources are present simultaneously in tfie INTEGRAL field-of- 
view. 

We used HEASoft Version 6.11 and selected tfie energy 
ranges 5-60 keV for PCA and 20-145 keV for HEXTE for 
tfie RXTE data. In tfie PCA spectra we excluded data below 
5 keV due to an instrumental background feature at tfie xenon 
L edge around ~ 4.8 keV, wfiicfi is poorly modelled and migfit 
introduce systematic eff'ects to our spectral fits (Rotfiscfiil d et al.l 
120061: ICabanero etsl\ 12008) We added a systematic error of 
1% for PCA and 1% for HEXTE data as recommended by tfie 
instruments teams. A systematic error for tfie HEXTE data is 
required due to tfie stop of tfie rocking motion of tfie HEXTE 
Cluster A in 2006 and HEXTE Cluster B in December 2009. 
Since tfien, HEXTE Cluster A is fixed on-source wfiile HEXTE 
Cluster B measures tfie background, w fiicfi is tfien converted to 
an adequate b ackground for Cluster A dPottscfimidt et alJl2006l) . 
Similarly to Rotfiscfiild et al.l (1201 ll) . we used tfie recorn model 
in XSPEC to account for uncertainties in tfie coiTection factor of 
tfie background. 

In order to select time intervals for tfie pfiase-resolved and 
pulse-to-pulse analyses, we extracted ligfit curves for eacfi ob- 
servation. We tfien converted tfie pfioton arrival times to tfie solar 
system barycenter and corrected tfiem for t fie orbital motion o f 
the binary system using the ephemeris from lFinger et al.l(l2006l) . 
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The pulse period was determined with the epoch-folding tech- 
nique. The pulse phase zero corresponds to the minimum in the 
pulse profile (see Fig.[8]i. 

3. Spectral analysis 

3.1. Phase averaged spectra 

The spectra of A 0535+26 for various luminosity states were 
best fitted by a power-law with a photon index F modified 
by an exponential roUofF with the folding energy /iroiioft (the 
cutof fpl model in XSPEQ: 

F{E) = const ■ E-^ ■ e--E/£r„iioif ^ 

We included an emission line at around 6.4 keV to model 
the Fe Ka fluorescence line. Additionally, the residuals of the fit 
require inclusion of a an absorption line around 40 keV (Fig. |2]l 
which we interpret as a CRSF (see Introduction). To account for 
this feature, we included a multiplicative absorption line with a 
Gaussian optical depth profile (the XSPEC gabs model) with a 
centroid energy of ~44.4keV: 



F(E) - exp 



' eye 



^/2no 



■ exp 



(£" £'cyc ) 



(1) 



where E^yc is the cyclotron line centroid energy, r^yc and cr^yc 
characterize the line strength and width respectively. The RXTE 
observations also indicate the presence of the first harmonic line 
around 104keV. The best-fit parameters for the INTEGRAL and 
RXTE observations are summarized in Table |2] Figure |3] shows 
the photon index F, the roUofF energy /SioUoff and the cyclotron 
line energy E^yc as a function of 10-lOOkeV flux in Crab units. 
The latter is obtained from a direct comparison of the source flux 
with that of the simultaneously observed Crab nebula. Because 
the first harmonic in the RXTE spectra is not statistically signif- 
icant and to be consistent with the INTEGRAL spectral model, 
only one Gaussian absorption line was included in our spectral 
model, accounting for the fundamental cyclotron line. 

The photon index decreases from ~ 1 .0 for low luminosities 
to ~0.5 for the outburst maximum, indicating a hardening of the 
power-law continuum with flux. The rolloff energy /SioUoif , how- 
ever, decreases from ~24 to ~16keV from the lowest to high- 
est luminosity. That is, the part of the continuum around /Si-oUoff 
softens with flux. No correlation of the cyclotron line energy 
with the source luminosity is observed. The cyclotron line en- 
ergy, however, varies irregularly between ~43 keV and ~46 ke V. 
To rule out any coupling between /SroUoft and F, we produced 
contour plots of those two parameters (Fig. 0. Since the evolu- 
tion of F indicates a hardening of the spectrum with increasing 
flux while the variation of /SioUoft suggests the opposite behav- 
ior, we calculated the hardness ratio to study the spectral vari- 
ations independently of the spectral model. Figure |5] shows a 
ratio of the 5-20 keV flux to the 20-50 keV flux as a function 
of the 10-100 keV flux. The observed decrease of the ratio im- 
plies an overall hardening of the spectrum with flux. A simi- 
lar behavior is reported in Caballero et al. (2013) for the nor- 
mal double-peaked outburst of August 2009. The results of this 
analysis, i.e. no correlation for the cyclotron line energy and a 
spectral hardening with flux, seem to support our results in the 
analyzed luminosity range. Following lAcciari et al ] (l2011h . we 
tried to model the spectra with the compTT model in XSPEC. We 
find that although the compTT model provides a good description 
of the X-ray spectrum below ~35 keV (the energy range used in 



lAcciari et al.l201 lb . it produces strong residuals when applied to 
the overall energy range evaluated in this paper, thus making the 
compTT overall fit formally unacceptable. Nevertheless, we ana- 
lyzed the evolution of the cyclotron line energy as a function of 
flux using compTT and found it to be consistent with the results, 
obtained with our phenomenological model. 

3.2. Pulse-phase resolved spectroscopy 

To select pulse-phase intervals, we produced pulse profiles 
whi ch are s imilar to those reported earlier (e.g. Caballero 200^ 
and ICaballero et al.l |201 Ibl) and do not show strong variation 
with flux. We then extracted spectra in 10 equal phase bins for 
each observation. That is, we performed pulse-phase resolved 
spectroscopy separately for different flux levels. Figure |6] sum- 
marizes the evolution of the photon index and the cyclotron line 
energy with pulse phase at different fluxes. For each pulse-phase 
bin separated by vertical dashed lines and marked by a num- 
ber from 1 to 10, left column, F and E^yc are shown as a func- 
tion of 10-100 keV flux in the middle and the right column, re- 
spectively. To increase the photon statistics in the observations 
at the lowest flux levels, we summed up the observations from 
the /A^TTiG/ML revolutions 0916 and 0917. For each pulse-phase 
bin, we performed a linear fit of the F-flux and /icye-Aux depen- 
dencies shown by the solid lines. The slopes of the linear fits 
with the corresponding Icr uncertainties are shown in the bot- 
tom panels of the left column. The value of the slope with the 
associated uncertainty gives a rough estimate of the significance 
of the linear correlation between the parameters. 

F and ficye show systematic variations with pulse phase and 
flux. The evolution of F as a function of flux, for individual 
pulse-phase intervals, indicates a hardening of the power-law 
part of the continuum with flux in almost all phase bins. In the 
pulse phases between 0.7 and 1.0 there is a strong indication 
(slope/error >3 cr) for a positive correlation of the cyclotron line 
energy E^yc with flux (bottom panel of left column in Fig.|6]l. To 
estimate quantitatively the significance of the non-zero slopes, 
we performed a F-test, which yields a value P of the probabil- 
ity that the improvement in x~ obtained by moving from a fit 
of the dependences by a constant to a linear fit with non-zero 
slope is due to statistical flucuations. Lower values of P, thus, 
indicate larger significance. In case of the F-flux dependence, 
the P-values for most phase bins are at the level of a few per 
cent indicating marginal significance for individual bins. For the 
cyclotron line energy dependence on flux in the phase interval 
0.7-1.0., i.e. the phase bins 8, 9, and 10, the P-values are 0.027, 
0.009, and 0.001, respectively. 

To check possible influence of any systematic differences 
in calibration between the INTEGRAL and RXTE instruments 
on our results, we re-calculated the significances using only the 
INTEGRAL data points. For the dependence of F on flux, the P- 
values remained at a level of a few per cent. For the E-finx de- 
pendence in the phase bins 8, 9, and 10, the P-values are found 
to be 0.08, 0. 17, and 0.05, respectively. The values for the slope, 
however, did not change considerably. We find, therefore, no in- 
dication that the measured correlations are driven by systematic 
differences between the instruments onboard the two satellites. 

As before, to check the stability of our results with respect 
to the choice of the continuum model, following Acci ari et al.l 
(1201 Ih we fitted the spectra using the thermal Comptonization 
model compTT in XSPEC. Even though the spectral fits are for- 
mally unacceptable (similar to the case of phase-averaged anal- 
ysis, see Section 3.1), the evolution of the cyclotron line energy 
with flux is well reproduced. 
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Table 2: Fit parameters for the INTEGRAL observations 0912, 0914, 0915, 0916 and 0917 and the KXTE observations 95347-02- 
01-00, 95347-02-01-01 and 95347-02-01-02. The errors are at 1 cr confidence level. 
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^rolloff 


p 








INTEGRAL 0912 
INTEGRAL 0914 
INTEGRAL 0915 
INTEGRAL 0916 
INTEGRAL 0917 


30+" " 

-0.13 

65+" " 

"•"-'-0.07 

n QQ+0.08 

"•^^-0.09 

1 or""** 

^•"'^-0.09 


17 91 +0-44 
^ '-^^-0.48 

15 32+"-" 

^-'•-^ -0.90 
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4^ qq+0.41 

^-'•''^-0.41 

44 31+0.70 

^^•-^^-0.67 

45.83!i;°« 
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^^•'-'-0.80 
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13 16+-"* 
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16.82^3-24 


1.1/110 
0.9/110 
0.9/110 
0.8/110 
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/fXr£ 95347-02-01-00 
RXTE 95347-02-01-01 
ifXre 95347-02-01-02 
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43+" °2 
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"■^'-0.02 


16 24+° " 

15 95+° " 
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44 94+0.29 

43 QO+0.29 

-0.29 

43 3r"-33 


QQj+0.36 

9 80+" 3' 

'■°"-0.34 

9 16+"^' 

'■^"-0.40 


8 69+° '2 

°-"'-0.49 

9 28+° ''' 
9 00+° '^3 


1.2/152 
1.5/152 
1.3/152 
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Fig. 2: The unfolded spectrum of A0535H-26 from the 
INTEGRAL revolution 0912 (a), the residuals of the fit 
with the model without the 6.4 keV Fe emission and the 
cyclotron absorption lines (b), the residuals after inclu- 
sion in the model of the cyclotron absorption line (c) and 
the Fe emission line (d). 



1.2 

1.0 

I 0.8 

0.6 

I 0.4 

& 0.2 

0.0 
26 

24 

¥22 

t 20 

I 18 



Luminosity in lO'' crg/s [10-100 keV| 
0.2 0.6 1.0 



i 



i 



16 - 
14 

47 

46 

45 

44 

43 

42, 



Fig. 3: Photon index F, rolloff energy /iioiioff and cy- 
clotron line energy E^yc as a function of the source flux 
in the 10-100 keV band in units of Crab and the lumi- 
nosity in erg/s for a distance of 2kpc (upper X-axis). 
The INTEGRAL observations are marked with circles, 
the RXTE observations - with triangles. The errors are 
at Icr confidence level. 
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Fig. 4: A contour plot of the folding energy Zsroiioff and 
the photon index F for all five INTEGRAL observations. 
The contour levels are at;if^j^ -H 1 .0 (corresponds to 68% 
confidence for one chosen parameters) and + 2.3 
(corresponds to 68% confidence for two chosen parame- 
ters). 
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Fig. 5: The 5-20 keV to 20-50 keV flux ratio as a func- 
tion of the source flux in units of Crab. The graph indi- 
cates a hardening the spectrum with flux. 
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Fig. 6: The best-fit parameters of the pulse-phase resolved spectra at different flux levels with 1 cr errors. Left column: In the top 
panel, the averaged and normalized pulse profile of the INTEGRAL observation during rev. 09 1 2 (3-80 ke V, black line) and the RXTE 
observation 95347-02-01-01 (3-60 keV, green line) are shown. The pulse-phase intervals for which the spectra were extracted are 
indicated by vertical dashed lines and numbered from "1" to "10". Within each bin, we analyzed the spectral parameters as a 
function of flux. The middle and bottom panels represent the slopes with associated Icr-uncertainties obtained from the linear fits 
to the F-flux and /icyc-flux dependencies in each pulse-phase bin which are shown in the middle and right columns. Middle and 
right columns: The photon index F and the cyclotron line energy E^yc as a function of the 10-100 keV flux for each of the ten phase 
intervals marked by the numbers. The solid lines show the linear fits to the data points. The RXTE observations are marked with 
triangles, while the INTEGRAL observations - with circles. 



In the following, the two peaks of the double-peaked pulse 
profile are referred to as the main and secondary pulse. To char- 
acterize spectral variations within larger pulse-phase intervals, 
namely in the main and secondary pulses, we extracted spectra 
for two broad pulse-phase intervals for each observation. The 
first interval, 0.1-0.5, contains the main peak, the second one, 
0.6-1 .0, - the secondary peak. The best-fit parameters, as a func- 
tion of flux, for the two phase intervals are shown in Fig.|7] To 
characterize the hardness of the continuum by the photon index 
F only, we fixed the rolloff energy /iroiiotf to its averaged value of 
16 keV for the upper plot of Fig. Q The spectrum is harder and 
the cyclotron line energy lower for the main pulse. We observe 
an indication for a spectral hardening and a positive trend of the 
cyclotron line energy with flux for the secondary peak. The val- 
ues for the slopes of the linear fit together with the results of the 
linear correlation analysis are shown in Table |3] 

3.3. Pulse-to-pulse analysis 

We performed a spectral analysis of the source on a pulse-to- 
pulse basis, similar to Klochkov et al. (2011). To identify indi- 
vidual pulsations, we extracted light curves for all I NTEG RAL 
and RXTE observations. Following Klochkov et al ] (120111) . we 
have chosen the secondary peak of the pulse profile (Fig. |8) for 
our analysis, to be able to compare our results. In the pulse pro- 
files above ~ 20 — 30 ke V the secondary peak is higher than the 
first one, thus, providing better statistics for the measurements 
of the cyclotron line around 50keV. For each individual pulsa- 



tion cycle, we calculated the mean ISGRI (for the INTEGRAL 
observations) and PCA (for the RXTE observations) count rate 
in the selected pulse-phase interval, which we refer to as height 
or amplitude of individual pulses. We then grouped pulses of a 
similar pulse height together to accumulate spectra as a function 
of pulse height. 

For the pulse-to-pulse analysis of the INTEGRAL observa- 
tions, we only used the data from the revolution 0912 taken 
during the maximum of the outburst, which provided the best 
photon statistics. Strong pulse-to-pulse variability can be seen 
clearly in Fig. |9] where the 20-120 keV ISGRI light curve, to- 
gether with the repeated averaged pulse profile and the selected 
pulse phase interval, is shown. The distribution of the measured 
amplitudes of individual pulses for this observation is shown in 
the left panel of Fig.[TOl We accumulated spectra in five pulse- 
height bins: 100-210, 210-235, 235-260, 260-290, and 290- 
450 ISGRI counts/s. 

For RXTE, we analyzed all three observations together as 
they were performed at a similar luminosity levels, close to the 
maximum of the source's outburst. The selected pulse height in- 
tervals for which the spectra have been accumulated are 850- 
1040, 1040-1210, 1210-1310, 1310-1480, and 1480-1850 PCA 
counts/s. The distribution of pulse heights for the RXTE data is 
shown in the right panel of Fig.fTOl 

The INTEGRAL and RXTE observations were analyzed sep- 
arately because of possible systematic effects of non-peifect 
cross-calibration between the two observatories. The best-fit 
photon index and cyclotron line energy as a function of pulse 
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Table 3: The slopes of the linear fits to the dependence of the spectral parameters on flux with the corresponding Pearson correlation 
coefficients and the associated two-sided null-hypotheses probabilities (a lower value indicates higher significance of the correla- 
tion). The uncertainties for the slopes are at 1 cr confidence level. The parameters from the pulse-phase resolved analysis of the phase 
intervals 0.1-0.5 (main peak) and 0.6-1.0 (secondary peak) are shown in the top section (for F, iiroiiotf was fixed to 16keV). In the 
bottom section, the values from the pulse-to-pulse analysis are shown. For the fit of F we fixed /iroiioff - 15.5 keV for INTEGRAL 
and /iioUoff = 15.2 keV for RXTE. For comparison, the values from Klochkov et al. (201 1) are also shown. 



Observation 



parameter linear fit 



correlation coefF. probability 



Pulse-phase-resolved analysis 

phase 0.1-0.5 

phase 0.6-1.0 



r 



-0.00 ± 0.03 (Crab unit)-' 
-0. 15 ± 0.70 keV (Crab unit)-' 
-0.12 ±0.03 (Crab unit)-' 
2.09 ± 0.95 keV (Crab unit)-' 



-0.28 
-0.57 
-0.83 
-0.60 



Pulse-to-pulse analysis 

INTEGRAL 

RXTE 

values from Klochkov et al. (2011) 

INTEGRAL 

RXTE 



(-5.79 ± 1.66) X 10-*/(ISGRI cts/s) 
(0.69 + 6.78) X lO--'' keV/(ISGRI cts/s) 
(-9.12 ± 1.84) X 10-V(PCA cts/s) 
(2.02 ± 1.26) X lO--'' keV/(PCA cts/s) 

(-8.71 ± 0.86) X 10-V(ISGRI cts/s) 
(1.36 ± 0.53) X 10-2 keV/(ISGRI cts/s) 
(-1.54 ± 0.19) X ^-^(PCA cts/s) 
(1.92 ± 0.92) X 10-3 ]ieV/(PCA cts/s) 



—0.96 
—0.17 
—0.99 
-0.78 

-0.96 
+0.94 
-0.98 
+0.70 



-0.51 
-0.14 
-0.01 
-0.12 



-0.01 
-0.78 

- 6.89 X lO-'' 
-0.12 

0.01 
0.02 

4 X lO-"* 
0.19 



height are shown in Fig.[TT] We fixed the rolloff" energy ZiroUoff to 
15.5 keV for the INTEGRAL and 15.2 keV, for the RXTE data 
for the plot of the photon index F. A significant decrease of 
the photon index (hardening of the power-law continuum) with 
pulse height is indicated by the data from both satellites. The 
cyclotron line energy, however, does not show any considerable 
evolution with pulse amplitude in the INTEGRAL data and an in- 
dication of a positive coiTelation with pulse height in the RXTE 
data. The results of the corresponding formal linear coiTelation 
analysis are shown in Table |3] Although the statistics of our data 
does not allow us to confirm the po sitive correlation be tween 
Ecyc and pulse amplitude reported in Klochkov et al.l (1201 1) (al- 
though, some support to the correlation is leaned by the RXTE 
data), the measured slopes are in agreement, with in uncertain- 
ties, with those reported bv Klo chkov et al.l (1201 ih and repeated 
in the bottom section of Table |3] 

4. Discussion 

We have analyzed the phase-averaged, phase-resolved and pulse- 
to-pulse spectra of A 0535+26 as a function of the source lumi- 
nosity. In all types of the analysis, the continuum gets harder 
with increasing flux. In some pulse-phase intervals, the data 
show indications of a positive Zicyc-flux correlation, although 
with marginal significance (~3 cr). The pulse-averaged centroid 
energy of the cyclotron line measured over longer time scales 
(days) shows some iiTegular variability but no clear dependence 
on flux, which is consistent with previous observations of this 
source. In the following, we discuss our results in more detail. 

4.1. Long-term spectral variations witli luminosity 

As mentioned in the introduction, a relation between the 
cyclotron line centroid energy and luminosity was reported 
for a number o f accr eting pulsars. Since the work of 
iBasko & SunvaevI (Il976h . it is known that the geometry of 



the emitting structure (column), on/above the polar caps of an 
accreting neutron star, must be different for the luminosities 
above and below a certain critical value L^ which depends on 
the pulsar's properties. In the following, accreting pulsars with 
luminosities below L^ are referred to as sub-critical sources 
while those with luminosities exceeding L^ - as super-critical. 
According to this picture, V 0332+53 and 4U01 15+63, which 
show a negative E^yc - L coiTelation and a softening of the 
spectrum with flux, belong to so-called super-critical sources, 
whose luminosities are above L^. In such sources, infalling mat- 
ter is most probably decelerated in a radiative shock, whose 
height is believed to increase with L, i.e., drift towards an area 
with a lower B-field strength. Her X-1 and GX 304-1, where a 
positive Ecyc - L correlation and a hardening of the spectrum 
with flux are observed, are the examples of sub-critical sources 
(L < Lc), where infalling matter is stopped by the Coulomb drag 
and collective plasma effects, rather than in a radiative shock. 
In such sources, the height of the emitting region decreases 
with increasing luminosity, owing to a corresponding incre ase 
in ram pressure of the infalling material (IStaubert et al.l20()7l and 
Becker et al. 2012). 

The newest theoretical calculations by iBecker et al.l (1201 2h 
predict for A 0535+26, a value of the critical luminosity L^ ~ 
6.78 X 10^^ erg/s, which is ~5-6 times larger than the peak flux 
of the outburst analyzed here (see the upper X-axis of Fig.|3). 
Thus, A 0535+26 must have been accreting in a sub-critical 
regime during the outburst. Though, the key indication of this 
regime, namely the positive correlation between the CRSF en- 
ergy and flux, is not firmly established in this source. Its X-ray 
continuum shows, however, a hardening with increasing flux, 
which is also observed in the other two sub-critical sources. Her 
X-1 and GX 304-1 (the super-critical sources V 0332+53 and 
4U 01 15+63 show an opposite behavior of the spectral contin- 
uum, see above). If the accretion regime in A 0535+26 is similar 
to one in Her X-1 and GX 304-1, the expected positive E^yc - L 
correlation must be smeared due to a special geometry of the 
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Pulse phase 

Fig. 8: The averaged pulse profile accumulated with 
INTEGRAL during its revolution 0912 (3-80 keV, continuous 
Une) and with RXTE in the observation 95347-02-01-01 (3- 
60keV, dotted line). The shaded area indicates the pulse-phase 
interval used for the pulse-to-pulse analysis. 
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Fig. 10: Distributions of individual pulse heights 
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Fig. 9: The 20-120 keV ISGRI Hght curve of the INTEGRAL 
observation during rev. 0912 (black data points) together with 
the overplotted repeated averaged pulse profile (solid curve). 
The selected pulse-phase interval is indicated by shaded ar- 
eas. The pulse-to-pulse variability can clearly be seen. The 
measured amplitudes of individual pulses (see text) are shown 
with magenta circles. 
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measured with INTEGRAL (left) and RXTE (right). 
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Fig. 11: Left: The results of the pulse-to-pulse analysis of the INTEGRAL observation of rev. 0912. The photon index F with fixed 
£^roiiofF=15.5keV and the cyclotron line energy ficyc are plotted as a function of the individual pulse height. The vertical error bars 
indicate Icr-uncertainties. The horizontal eiTor bars correspond to the pulse-height bins used to accumulate the spectra. The soli3 
lines indicate linear fits to the data. Right: The same but for the RXTE observations. In the data for the upper panel /iroUoff was fixes 
to 15.2 keV. 
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Fig. 7: Best-fit spectral parameters as a function of flux for 
the main peak (phase interval 0.1-0.5) and the secondary peak 
(phase interval 0.6-1.0) of the pulse profile. The data points and 
the linear fit (solid line) corresponding to the main peak are 
shown in black. The data corresponding to the secondary peak 
are shown in magenta. As in previous figures, the INTEGRAL 
observations are marked with circles, the RXTE observations - 
with triangles. For the upper plot, /SroUoff was fixed to 16.0keV 
(see text). For the bottom plot all fit parameters are kept free. 
The corresponding values for the linear fit and the significance 
are shown in Table [3] 



emitting structure and/or a special orientation of the system with 
respect to the observer's line of sight. Alternatively, A 0535+26 
might represent yet another, "third" accretion regime which 
yields very little, if any, variations in the height of the emit- 
ting region above the star surface (i.e., in the measured cyclotron 
line energy) in the observed range of luminosities. Indeed, while 
the luminosities of Her X-1 and GX 304-1 are below, but rel- 
atively close to th e corresponding values of (see Fig. 2 in 
lBeckeret al.120121) . A 0535+26 is far below Lc- Thus, the config- 
uration of the accretion column in this source might differ sub- 
stantially from those in the other two sub-critical pulsars. These 
considerations seem to be supported by the lat est observations 
of a g iant outburst of A 0535+26, analyzed bv ICaballero et al.l 
(I2012h . when the source luminosity approached Lc (i.e., a lumi- 
nosity range of Her X- 1 and GX 304- 1 ). In this preliminary anal- 
ysis, a few measurements of iScyc, corresponding to the highest 
luminosities, indeed show an in dication of an increase towards 
higher fluxes (Fig. 1 of Caballer o et al.ll20T2b . Further observa- 
tions are, however, necessary to confirm this effect. 

4.2. Pulse-phase dependence of the spectrum 

It is necessary to reconstruct the geometry of the rotating neu- 
tron star and its emitting regions to understand the results of the 
phase-resolved spectral analysis. Such a reconstruction can be 



made based on the observed pul se profiles of a pulsa r using, e.g., 
the decomposition method by Kraus et aP (11995). Thi s tech- 
nique h as recently been applied to A 0535+26 by Caballer o et al.l 
(|201 Ibl) . As a result, the pulse profiles of the source were repro- 
duced in a model, including a hollow column, a thermal halo 
and taking into account relativistic light deflection in the strong 
gravitational field of the neutron star Derived from the results of 
Caballero et al. (2011b), the observed pulse profile can be pro- 
duced by one magnetic pole passing through the observer's line 
of sight, whereas the second pole is never seen directly from the 
top but only through the gravitationally bent radiation from the 
lateral emission of the accretion column and its thermal scatter- 
ing halo. 

Since we o b served similar pulse profiles as those used by 
ICaballero et al.l (|201 Ibl) . we adopt the reconstructed geometry 
for our data. In our analysis, we observe a slightly harder spec- 
trum for the main pulse with respect to the secondary pulse. 
Magnetized plasma gets optically thin for photons (E< E^yc) 
which propagate parallel to the magnetic field lines. Thus, the 
main peak is probably formed by harder photons escaping the 
accretion column from top, along the field lines. About half a 
phase later, we probably see the lateral emission of both the ac- 
cretion column and their thermal halos. 

In our phase-resolved analysis, we observe variations of the 
cyclotron line energy with pulse phase in range ~42-50keV, 
which is probably caused by the changing viewing conditions of 
the scattering region and a dependence of the magnetic scatter- 
ing cross sections on the angle between the m agnetic field and 
the line of sight dHarding & Daughertvl[T99TI) . Variations over 
pulse phase are also seen in the slope of the linear fits to the 
Ecyc and F values as a function of flux (Fig.|6ll. While a signif- 
icant decrease of F with flux is present at most pulse phases, a 
considerable (positive) ZScyc-fux correlation is only seen during 
the secondary peak (Figs.|6] and |7]i. According to the geomet- 
rical picture described above, the emerging radiation is emitted 
roughly perpendicular to the axis of the accretion column(s), that 
is, perpendicula r to the field lines, in the secondary peak. As 
discussed in e.g.'Ferrig no et al.l (1201 ll) . at such angle, the reso- 
nant scattering features are expected to be deeper and narrower 
Therefore, in the secondary peak, the viewing conditions of the 
line-forming region might be favourable for measuring E^yc vari- 
ations with flux. In the pulse-averaged spectral analysis, the de- 
pendence is smeared out. 



4.3. Pulse-to-pulse variability 



As we noted in Sect. l3.3l the results of our pulse-to-pulse 
analysis, namely the spectral variations with the amplitude of 
individual puls e s, are in agreement with those reported by 
iKlochkov et alJ ll2011h on the basis of the observations of the 
source's outburst in 2005 with better photon statistics. With our 
data, we have reproduced the hardening of the power-law with 
pulse amplitude. An indication of a positive correlation of the 
cyclotron line energy with flux can only be seen in the RXTE 
data (Fig.fTTI). If the positive correlation between the line energy 
and flux on the pulse-to-pulse time scale in A 0535+26 is real, 
the question remains: why the correlation is not observed on the 
pulse-averaged spectra accumulated over longer time scales? A 
qualitative explanation of this discrepancy might be the follow- 
ing. Strong pulse-to-pulse variability of the source's flux leads to 
the relatively broad distribution of the amplitudes of individual 
pulses - more than a factor of two, as indicated in Fig.[TO] Thus, 
pulse-averaged spectra accumulated over a longer time include 
contributions by pulses of very different heights. The pulse-to- 
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pulse variability of the source's spectrum must lead to a substan- 
tial smearing effect, especially in case of a line-like feature. If 
A 0535+26 indeed accrets in a regime with very little (if any) 
variation in the height of the emitting region with flux, as dis- 
cussed in Sect. l4.1l the reduced sensitivity of the pulse-averaged 
spectroscopy to the /Scyc-flux coiTelation might indeed lead to 
the observed discrepancy between the pulse-averaged and pulse- 
to-pulse results. 

5. Summary and conclusions 

In this work, we present the spectral analysis of the observa- 
tional data taken with RXTE and INTEGRAL on the Be/X-ray 
binary system A 0535+26, during its outburst in April 2010. We 
studied the source's broad-band X-ray spectrum as a function 
of luminosity on the pulse-averaged, pulse-phase resolved, and 
pulse-to-pulse basis. In all cases, a hardening of the X-ray con- 
tinuum with flux is observed. This behavior is similar to that of 
other sub-critical sources. Her X-1 and GX 304-1. An indica- 
tion of a positive correlation between the cyclotron line energy 
and flux (a characteristic feature of a sub-critical source) is only 
found in some pulse-phase resolved and pulse-to-pulse spectra 
(Fig.|6]andTable[3]l. 

Based on the spectrum-flux dependence of A 0535+26, 
which is similar to that of other sub-critical sources, and on the 
source's luminosity, which is well below the critical luminosity 
Lc, we conclude that the pulsar operates in a sub-critical accre- 
tion regime, where the radiative shock, if formed, is not capable 
of decelerating the accreted plasma to rest. Instead, the in-falling 
matter is most probably stopped by the Coulomb drag and col- 
lective plasma effects. A similar behavior was reported in ear- 
lier wo rks on this source, see e.g. Caballero et a l. (2007, 20()8, 
|20TT^ . 

The lack of a firm detection of a Ecyc-Aux correlation, which 
is expected in this accretion regime, we tentatively attribute to 
a substantially lower ratio of the source luminosity to the crit- 
ical one, compared to other sub-critical sources. At such a low 
luminosity level, the amplitude of variation in the height of the 
line-forming region with luminosity (mass accretion rate) might 
be substantially narrower than that in Her X-1 and GX 304-1. 
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